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We study the perturbation modes of rotating superfluid ellipsoidal figures of equilibrium in 
the framework of the two-fluid superfluid hydrodynamics and Newtonian gravity. Our calculations 
focus on linear perturbations of background equilibria in which the two fluids move together, the 
total density is uniform, and the densities of the two components are proportional to one another, 
with ratios that are independent of position. The motions of the two fluids are coupled by mutual 
friction, as formulated by Khalatnikov. We show that there are two general classes of modes for 
small perturbations: one class in which the two fluids move together and the other in which there is 
relative motion between them. The former are identical to the modes found for a single fluid, except 
that the rate of viscous dissipation, when computed in the secular (or "low Reynolds number") 
approximation under the assumption of a constant kinematic viscosity, is diminished by a factor /jv , 
the fraction of the total mass in the normal fluid. The relative modes are completely new, and are 
studied in detail for a range of values for the phenomenological mutual friction coefficients, relative 
densities of the superfluid and normal components, and, for Roche ellipsoids, binary mass ratios. We 
find that there are no new secular instabilities connected with the relative motions of the two fluid 
components. Moreover, although the new modes are subject to viscous dissipation (a consequence 
of viscosity of the normal matter) , they do not emit gravitational radiation at all. 



I. INTRODUCTION 



The problem of the equilibrium and stability of rotating neutron stars is encountered in various astrophysical 
contexts, ranging from the limiting frequencies of rapidly rotating isolated millisecond pulsars, emission of gravitational 
waves in neutron star-neutron star or neutron star-black hole binaries, to the generation of 7-rays in the neutron star 
mergers and X-rays in accreting systems Considerable current interest is attached to the problem of neutron 
star binary inspiral, which would be the primary source of gravitational wave radiation for detection by future laser 
interferometers. Such a detection, apart from testing the general theory of relativity, potentially could provide useful 
information on the equation of state of superdense matter. Also the stability criteria for rapidly rotating neutron 
stars are essential for placing firm upper limits on the frequencies to which millisecond pulsars can be spun up thereby 
constraining the range of admissible equations of states. 

High precision, fully relativistic treatments of rapidly rotating isolated neutron stars and binaries comprising two 
neutron stars have become available in recent years Nevertheless, the development of simpler models that provide 
a fast and transparent insight into the underlying physics is needed when the basic set of equations is modified to 
include new effects. 

A systematic framework for the treatment of the equilibrium and stability of rotating liquid masses bound by self- 
gravitation in the Newtonian theory is contained in Chandrasekhar's Ellipsoidal Figures of Equilibrium (hereafter 
EFE). The tensor virial method, developed most extensively by Chandrasekhar and co-workers, transforms the local 
hydrodynamical equations into global virial equations that contain the full information on the structure and stability 
of the Newtonian self-gravitating system as a whole. The method describes, in a coherent manner, the properties of 
solitary ellipsoids with and without intrinsic spin, and ellipsoids in binaries subject to Newtonian tidal fields. It is 
especially useful for studying divergence-free displacements of uniform ellipsoids from equilibrium, in which case the 
each perturbed virial equation yields (in the absence of viscous dissipation) a different set of normal modes. 

Recent alternative formulations of the theory of ellipsoids are based mainly on either the energy variation method 
|^,|| or the afhne star model or the (Eulerian) two potential formalism Q. A large class of incompressible and 
compressible ellipsoidal models has been studied using these methods [^-0 . The energy variation method employs the 
observation that an equilibrium configuration is possible if the energy of the ellipsoid is an extremum for variations 



* Postal address. 



1 



of the ellipsoidal semiaxis at a constant volume; the ellipsoidal figure is stable only if the energy is a true local 
minimum. In the affine star method the figures are described by a time-dependent Lagrangian as a function of a 
deformation matrix and its derivative. The structure of the star at any particular time is related conformally to the 
initial unperturbed sphere via a quadratic form constructed from the deformation matrix. 

The purpose of this paper is to extend previous studies to a treatment of the oscillation modes of ellipsoidal figures 
of equilibrium that contain a mixture of normal fluid and superfluid. Many-body studies of the pair correlations in 
neutron star matter show that the baryon fluids in their ground state form superfluid condensates in the bulk of 
the star. The superfluid phases, in the hydrodynamic limit, can be treated as a mixture of superfluid condensate 
and normal matter. The superfluid rotation is supported by the Feynman-Onsager vortex lattice state, and on the 
average leads to quasi-rigid body rotation of the superfluid component. The corresponding time-dependent two-fluid 
hydrodynamic equations are completely specified by two independent velocities for the superfluid and the normal 
component and corresponding densities of the constituents. The two fluids are coupled to one another by mutual 
friction forces, which we model phenomenologically according the prescription given by Khalatnikov p7[ |. Because 
of the additional degrees of freedom in this system, there are twice as many modes as for a single fluid. A natural 
question is whether the new modes affect the stability criteria previously deduced from studying perturbations of a 
single self-gravitating fluid. 

In our treatment of perturbations of a mixture of normal fluid and superfluid, we shall follow most closely Chan- 
drasekhar's formulation. However, since the basic equations of motion for the two fluids will include mutual friction 
between them, the system we study is inherently dissipative. Nevertheless, since the frictional forces only depend lin- 
early on the relative velocity between the two fluids and vanish in the background, where the two fluids move together, 
one can still derive relations that resemble Chandrasekhar's tensor virial equations. Because these equations include 
dissipation, we shall prefer to regard them as moments of the fluid equations, rather than tensor virial equations. In 
fact, we shall relegate the derivation of perturbation equations from these moment equations to the Appendix of the 
paper, and instead derive the necessary equations for the fluid displacements directly by taking moments of linearized 
equations of motion for the two fluid components. In this paper, we concentrate on two- fluid variants of the classical 
Maclaurin, Jacobi and Roche ellipsoids. 

Although our main aim is to access the oscillation modes and instabilities of neutron stars within the ellipsoidal 
approximation, the results obtained here may be of significance in other contexts (Ref. the epilogue). One example 
is the understanding of rapidly rotating nuclei in the spirit of the Bohr- Wheeler model of a charged incompressible 
liquid droplet In this case, the stability is determined by the competition between the attractive surface tension, 
the repulsive Coulomb potential, and the centrifugal stretching due to the rotation Another example is the 

stability of rotating superfluid liquid drops of Bose condensed atomic gases, where the stability is determined through 
an interplay among the pressure of the condensate, the confining potential of the magnetic trap and the centrifugal 
potential [jlTll . 

Previous work on the oscillations of superfiuid neutron stars concentrated mainly on perturbations of non-rotating 
or slowly rotating isolated neutron stars jl^-p^, and used methods entirely different from the one adopted here. The 
propagation of acoustic waves in neutron star interiors, including those related to the relative motion of neutron-proton 
superfluids, was studied by Epstein | [T^ , who found the compressional and shear modes related to short- wavelength 
oscillations of neutron star matter. The small-amplitude pulsation modes of superfluid neutron stars were derived 
by Lindblom and Mendell JTsf , who found that the lowest frequency modes were almost indistinguishable from the 
normal modes of a single fluid star. Their analytical solutions also reveal the existence of a spectrum of modes which 
are absent in a single fluid star. Subsequent work concentrated on numerical solutions for the radial and non-radial 
pulsations of the two-fluid stars and identified distinct superfluid modes in the absence of rotation . The r modes 
of slowly rotating two-fluid neutron stars have been derived by Lindblom and Mendell |l3|, who find that they are 
identical to their ordinary-fluid counterparts to the lowest order in their small-angular-velocity expansion. The linear 
oscillations of general relativistic stars composed of two non-interacting fluids in a non-rotating static background 
have been studied by Comer et al. . 

Our calculations allow arbitrary fast rotation, in the context of (incompressible) Newtonian fluid models. One may 
anticipate that the effects of superfluidity on oscillation modes, if any, should be affected by the underlying vortex 
structure of the rotating superfluid. In our treatment, dissipation arises because of the drag forces experienced by 
the vortex lines as they move through the normal fluid (and there is no dissipation if the drag force is zero). We 
ignore the motions related to the isospin degrees of freedom in the core of a neutron star and, hence, the mutual 
entrainment of the neutron and proton condensates, as well as forces arising from deviations from /3 equilibrium. The 
two-fluid equations used in the remainder of this work can adequately describe the mutual friction of a two-condensate 
fiuid in the core of a neutron star, since the entrainment effect renormalizes the effective superfiuid densities and the 
frictional coefficients, i.e. the phenomenological input in the two-fluid equations. While the ellipsoidal approximation 
to a superfluid neutron star is restrictive, it allows us to study the effects of vorticity on the oscillation modes of 
a self-gravitating star in a transparent manner, avoiding complications due to the star's inhomogeneity (multi-layer 
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composition) [l^ . 



II. PERTURBATION EQUATIONS 

The equations of motion for a mixture of two fluids may be summarized simply as 

dpc, d(t) 1 a|J^Xx|2 ^ „ ^ 

where the subscript a G {S, N} identifies the fluid component, and Latin subscripts denote coordinate directions; pa, 
Pa, and are the density, pressure, and velocity of fluid a, (j) is the gravitational potential, and F^^ is the mutual 
friction force on fluid a due to fluid (3. These equations have been written in a frame rotating with angular velocity 
Jl relative to some inertial coordinate reference system. The total time derivative operator 
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The gravitational potential (j) is derived from 

V^q} = V2(</)S + M = 47rG[ps(x) + pjv(x)]; (3) 

the individual fluid potentials (f>a obey 'S/'^(f>a = 47rGpQ. The two fluids are coupled to one another via the frictional 
force Faf3 which is antisymmetric on interchange of a and f3. For a normal-superfluid mixture 

FsN = -Fns = psws{/3'i'X(us - un)+ I3ux[ux{us - ujv)] - /3'V-(us - un)}, (4) 

where (3, (3' and (3" are coupling coefficients, and uis = ^"^5 = Vxu^; in components we have 

FsN,i = -Ps(^sf3ijiusj - UN,j), (5) 

where, from Eq. (Q), 

(3^j = /3% + P'e^j^iy,n + {13" - I3)v,vj. (6) 

The net rate at which this force does work is 

Us-Fsjv + Uat-Fats = -psLOs{PWx{\ls - UAr)P - f3"[i'-{\is - ujv)]^}; (7) 

there is no dissipation associated with the term proportional to j3' in Fq,^. Throughout this paper, we assume that 
(3, j3' and j3" are independent of position in the fluid mixture. In Eq. (||), we have neglected the effects of the 
vortex tension, and expressed the mutual friction force in terms of the phenomenological coefficients /3, (3' and (3" . 
While these parameters determine the macroscopic behavior of the fluid system, they are not the optimal ones for 
connecting microscopic parameters of the mixture to its macroscopic motion. Instead, the macroscopic results can be 
parametrized in terms of frictional coefficients ry and 77', which connect [3 and /?' to the drag on individual superfluid 
vortices via the relations 

^ ^ -qpsi^s /3' = 1 - P^^^^P'^^'^ ~ (8) 

772 + {psujs - rj'Y 772 + {psujs - v'Y ' 

The physical meaning of 77's is apparent from the equation of motion of a single vortex line 

PSUJS [(u5 - ul) X I/] - 7/ (ul - Un) + ri' [(ul - uw) x i^] = 0, (9) 

where is the vortex velocity. Equation (^) states that the Magnus force, which represents a lifting force due to the 
superflow imposed on the vortex circulation, is balanced by the viscous friction forces along the vortex motion (the 
term oc rf) and perpendicular to the vortex motion (the term <x rj'); these latter forces arise from the scattering of the 
normal quasiparticles off the vortex line[]. The characteristic dynamical relaxation time scale related to the vortex 
motion can be deflned as 



^The inertial mass of the vortex is neglected in the standard formulation of the two- fluid superfluid hydrodynamics |l7 
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r.^j^(^ + ^), (10) 

where (ws) is the superfluid circulation averaged over macroscopic scales; e.g. for uniformly rotating superfluid (ws) = 
2ri5, where is the superfluid rotation frequency. For fixed density p5, td — > oo asymptotically, when r\ ^ Ps^s 
(strong coupling limit) and psi^Js V (weak coupling limit). Its minimal value is attained when rj/psujs — 1- 

The relations (^ do not determine /?", which, if nonzero, implies friction along vortex lines, which would be possible 
if vortices oscillate or are deformed in the plane perpendicular to the rotation axis. Generally, we assume in this paper 
that (3" ^ (3 and (3' , but occasionally we retain nonzero (and not necessarily negligible) (3" to examine its effects on 
the modes. Moreover, if we assume that the drag force on a vortex line is principally opposite to its velocity relative 
to the normal fluid, then rj' <C 77, and the relationships between f3 and (3' and the drag force on vortex lines simplifies 
accordingly. 

Ultimately, we want to compute perturbations of Eq. (^ around some presumed background state. In general, 
the two fluids need not occupy the same volume, and we shall suppose that fluid a occupies a volume V^. However, 
we shall restrict ourselves to background states in which the two fluids occupy identical volumes and have densities 
Pq,(x) = /q,/o(x), where p(x) is the total density and does not depend on x. The perturbation equations derived 
below can be applied to nonuniform background states satisfying these conditions, but we shall only consider uniformly 
dense backgrounds, as was done by Chandrasekhar. Note, though, that we shall not assume that the two fluids must 
occupy the same volumes in the perturbed state. In fact, at least for low order perturbations, there are no nontrivial 
perturbations that leave the volumes of the two fluids identical. 

Following the example set by Chandrasekhar, we could take moments of the fluid equations (|^) to obtain tensor 
virial theorems of various orders, and then perturb them to find linear modes for uniform ellipsoids. We have derived 
the necessary moment equations in this way (see the Appendix), by analogy to Chandrasekhar's treatment for a single 
fiuid, but here we present a somewhat different (and possibly more transparent approach) to their derivation. We 
begin with the equation of motion for the displacement of fluid a from equilibrium, a direct generalization of Eq. 
(107) in Chap. 2, Sec. 14 in EFE: 



Pa ,2 — on o P'^ a ~^ Pa 

oxi oxi oxi oxi oxi 



2 V dxi 



9 O 

^Pa^ilm^^m ,, 

dt 



Po'l^- +Fcf3,^, (11) 



dxi \ dxi " dxi 2 dxi 

where the Lagrangian variation AaQ denotes the change in Q seen by a moving element of fluid a; in particular, it 
is easy to show that Ac0(x) separates into an Eulerian part, S(j), plus Ca.id4'/9xi, where 

S^^-GJ: //Vp,(x')e.,(x')|,(^) (12) 

is the Eulerian potential perturbation. We assume that the mutual friction force is Fap^i — in the background 
solution, i.e., either the two fluids are stationary in the rotating frame for the background or have identical fluid 
velocities in this frame. Our strategy will be to take moments of Eq. ( |ll|) by multiplying by appropriate factors of 
Xi and integrating over the (common) volume of the unperturbed background configuration. Note that this does not 
restrict the volumes of the perturbed fiuids in any way. 

Although the method we shall use to derive the perturbations resembles Chandrasekhar's, there is an important 
distinction due to the mutual friction force. Chandrasekhar's method of solution yields exact modes only in the 
dissipationless limit, where the modes of uniform ellipsoids are, respectively, linear, quadratic, cubic, etc. functions of 
the coordinates. Viscous terms would prevent exact solution in this manner, and one resorts to an approximation in 
which they are evaluated by substitution of the inviscid eigenfunctions (EFE, Chap. 5, §37(b)). However, it is possible 
to employ the moment method to find exact modes for the normal fiuid-superfiuid mixture coupled by mutual friction 
because Fap is a linear function of the velocity difference between the two fluids, not of their spatial derivatives 
(as is the case for viscous dissipation). By taking moments, we can derive the analogue of tensor virial theorems of 
various orders, but because the equations involve manifestly dissipative mutual friction forces, we prefer to think of 
these merely as moments of the original equations of motion. Of course, if we also choose to include the effects of 
viscous dissipation in the normal fluid, which we shall not do in this paper, we must resort to a low Reynolds number 



approximation, as was done in EFE. We discuss this briefly in Sec. IV 



Using the above result for the perturbed gravitational potential, we can rewrite Eq. (nil) as 
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" dt'^ dxi dxi ""'^ dxidxi dxi dxi 

— Pa^a,l^ ^ 1- Pa^a,l{^^Sil ~ ^i^l) + 2patilrn^rn ^"'^ + Fafi^i- (13) 



Then if we define 



we find 



oxiOxi at 



_ d£,+A dp ^ d^p dSp dS<j) 
dfi dxi dxi ^'^ dxidxi dxi ^ dxi 



0X^0X1 QAj 



where 5p = Sps + 5pN, and 



p d^^-,t _ dC-,i dp ^ ^ d^p 1 dSps ^ 1 dSpN 
dt^ dxi dxi ' dxidxi fs dxi /n dxi 

oxioxi dt 
-pu:s(l + ^\p.u'^. (16) 



In)' dt ■ 

Equation (|l^) is identical to what is found for a single fluid, and therefore contains the well-known modes documented 
by Chandrasekhar: if we define 

Vi-j = / d^xp£,+^,Xj, (17) 
Jv 

then we find 

- TTGp(2B,jV^j ~ a\h, ^a^u) ^ (18) 
^ 1=1 ^ 

where (511 = Jllg + (511 jv and all other quantities are defined exactly as in EFE. All of the new modes of a mixture 
of normal fluid and superfiuid are contained in Eq. dl^) for their relative displacements. One noteworthy feature of 



Eq. (16) is that the Eulerian gravitational potential does not appear. Consequently, the new normal modes of the 
system only depend on the unperturbed gravitational potential; for perturbations of homogeneous ellipsoids, only the 
coefficients Ai defined by EFE, Chap. 3, Eqs. (18) and (40), will appear. 

For the most part, we shall be interested in displacements that are linear functions of xi in this paper. For the 
homogeneous ellipsoids, we can find the new modes that result from the differential displacements of normal fluid and 
superfluid by taking the first moment of Eq. (p^), e.g. by multiplying by Xj and integrating over the unperturbed 
volume .0 If we define 

U,,j = / d^xp^^^.Xj, (19) 

JV 

then we find 



^In the nomenclature of EFE, this is the second order virial equation corresponding to Eq. (p^. 
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cPUi-j dUij 2 o o 

— = 2eilmilm 77^ h " Uij — iliilkUk 



dt 



1 ^ ^ij 



N 



V fs 



IN 



2^^GpAiU^J -ivsil 



fs\^ dUk 



f. 



N 



dt ' 



(20) 



where 



Slla = / d^xSpa- 



(21) 



To obtain Eq. (|2(]|), various surface terms can be eliminated using the conditions that Pa and AaPa = 5pa+£,a.idpa/ dxi 
vanish on the boundary; also, the equation of hydrostatic equilibrium for the unperturbed configuration, 







dp dcf) 



dxi 



^ dxi ^ dxi 



(22) 



must be used. It is straightforward to compute higher moments of Eq. (16). For example, taking its second moment^ 
by multiplying by XjX^ and integrating over the unperturbed volume gives 



^ i\jk 



dt^ 



6n 



Oik 



5n 



N,j 



V fs 



V fs fn 
2TrGpA,Uijk + {^"^Su - n,ni)Uijk 

dUi-jk 



fN 



dt 



where, by analogy to definitions in EFE for a single fluid. 



(23) 



i-jk 



d xpS^^^iXjXk 



Uijk = Ui:jk + Uj:ki + Uk;i 



SUa^k 



d XXkSpa- 



(24) 



Equations (EOf) and (03) only apply to unperturbed states that are static in the rotating frame. 



III. NEW MODES OF A TWO-FLUID MIXTURE 



In this section, we derive the characteristic equations for the normal modes for relative fluid displacements implied 



by Eq. (20) assuming Ui-j oc exp(Ai). We list the results separately for perturbations of Maclaurin, Jacobi and Roche 
ellipsoids. Note that the modes described by Eq. ( p^ ) are identical to those treated in EFE. 

The superfluid state of the matter does not affect the equilibrium figure.^ In equilibrium, Eq. ( pO| ) is satisfied 
trivially, for all virials Uij = in the absence of relative motion between the superfluid and normal components. The 
equilibrium figure follows from Eq. ( p^ by dropping the temporal variations and is identical to its non-superfluid 
counterpart. 



For irrotational ellipsoids we can set ujs — 2f2 in Eq. ( |20| ) and using the compact notations 

—? 7 : Pik — + T 

is JN J \ JN 



<5n=(^-^Y A. = fl + ^)A., (25) 



^This is the third order virial equation in the nomenclature of EFE. 
Here, we have taken pa — faP in the background state, which is a mathematically convenient idealization. More realistically, 
contributions from the pair condensation energy and the energy density of the superfluid vortex lattice, which distinguish the 
superfluid from the normal fluid, could play a role in both the equilibria and the perturbations. We need not restrict our 
attention to displacements for which the pressure perturbations of the two fluids are still proportional to one another, although, 
as argued in i: IV , for the solenoidal displacements considered here, this turns out to be the case. 
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we rewrite the Eq. (g6[) as 

- 2eara^-m ^ + ^^U,, - n,nkUkj + SU - 2^GpA,U,, - 2np,k • (26) 
For time-dependent Lagrangian displacements of the form 

U^^,t)^Ux,)e^\ (27) 

Eq. ( p^ ) becomes 

X^U,,j - 2eu,n^mXUi,j = +n^U^j - n,nkUkj + {ttpG)-^6,,6U - 2A,U^j - 2nxp,kUk;j; (28) 



here all frequencies are measured in the units (vrpG)^/^. Equation (^8|) contains all the second harmonic modes of the 
relative oscillation of Maclaurin and Jacobi ellipsoids, and only requires a slight modification for the application to 
Roche ellipsoids. 



A. Superfluid Maclaurin spheroid 

Next, we specialize Eq. ( p^ ) to Maclaurin spheroids, the equilibrium figures of a self-gravitating fluid with two 
equal semi-major axis, say oi and 02, uniformly rotating about the third semi-major axis (i.e. the ^3 axis). The 
sequence of quasi-equilibrium figures can be parametrized by the eccentricity = 1 — a|/af , with (squared) angular 
velocity fl'^ = 2€'^Bi3, in units of {wpGy^^. 

Surface deformations related to various modes can be classified by corresponding terms of the expansion in surface 
harmonics labeled by the indexes m. Second order harmonic deformations correspond to I = 2 with five distinct 
values of m, —2 < m < 2. The 18 equations represented by Eq. ( p6| ) separate into two independent subsets which are 
odd and even with respect to index 3. The corresponding oscillation modes can be treated separately. 



1. Relative transverse shear modes 



These modes correspond to surface deformations with |m| = 1 and represent relative shearing of the northern and 
southern hemispheres of the ellipsoid. The components of Eq. (Eg), which are odd in index 3, are 



A'f/3;1 



-2A3C/3i-2r!/3"AC/3.i, 

^3;2--2A3C/32-2r!/3"AC/3;2, 

x'^Ui,3 - 2nxu2;3 = {-2Ai + n^) U13 - 2iipxui,3 - 2np'xu2;3, 

X^U2,3 + 2nXUu3 = {-2Ai + f72) U23 - 2n(3XU2:3 + 2np'xui.,3. 

Note that because of the degeneracy in indexes 1 and 2 for the Maclaurin spheroid Ai = At. We sum Eqs. (p9| 
and (|30|), (^), respectively, and use the symmetry properties of Uij combined with Eqs. (|30|), (^). We find 

Ui3 



X^ + 2np"X )[X^ + 2VLPX + 2 A-' + 2Q.I3X M3 + A^ + 2Vl(3"X {2Ai - O^) 



-2nX{l-p') [X^ 



2n(3"X + 2A3 U: 



0, 



(a^ + 217/3" a) (a^ -f- 2nf3X^ + 2 (a^ + 217/3a) A3 + (^A^ + 20/3" A^ {2Ai - fl^) 
+2nXil - P') (X^ + 20/3"A + 2A3) Ui3 = 0. 



23 



It is instructive to consider first the limit of zero mutual friction, in which case Eqs. (p3|)-(|34|) reduce to 

A [A2 + {2Ai + 2A3 - n^)] Ui3 - 20 (A2 + 2A3) U23 - 0, 
A [a2 + {2Ai + 2A3 - 02)] U23 + 20 (a2 + 2A3) Ui3 = 0, 



(29) 
(30) 
(31) 
(32) 



(33) 



(34) 



(35) 



excluding the trivial mode A = 0. The characteristic equation can be factorized by substituting A = ia 
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The purely rotational mode cr = decouples only in the spherical symmetric limit where Ai 
symmetry is imposed then the characteristic equation is third order: 



(36) 

A3. If only axial 



cr^ ± 2r2cr2 ^ {Ai + A3) + f^^] cr T 4 A3 = 0. (37) 
Along the entire sequence parametrized in terms of the eccentricity the three modes are real N and are given by 



CTl 



2n , 



2VL 



i(s+ + s_) ± ^(s+ - s_), 



(38) 



where 



2r2(Ai - 2A3) - 



± 



27 3 

41^2 2(Ai+A3)-Sl2 



2A3 



0= 



27 



1/2 



(39) 



Three complementary modes follow from Eqs.(p8|)-(p9|) via the replacement — fi. The frictionless modes are real 
and are shown in Fig. 1. The two high frequency, frictionless modes are roughly twice as large as the transverse shear 
modes of ordinary Maclaurin spheroids. The third low frequency mode corresponds to the nearly rotational mode and 
indeed coincides with in the limits e — > and Ai — > A3, but not generally. In the dissipative case the characteristic 
equation is of sixth order: 



+ 



4A3 + 2f7(l - 13'Y + An^iP + I3"Y - 2(2Ai - f7^) 
IQA-ijirt + 8A3/?"0 + 8/3"(l - 13' fn'^ + lQi3^j3"n^ 

+16/3/3"2rj^ - 4/3rj(2Ai - VY) - 8/3"17(2Ai - Q^)] A^ 
AAl + 8A3l7(l - p'f + l&Aip^9? + iiAip^''^? + - p' f^l^ + 

-4A3(2Ai - O^) - 16/3/3" rj2(2Ai - - si3"^if{2Ai - if) + (2Ai - Vi^f 
l&Alpn + 16A3^"(1 - p'fif + 2,2A-iP'^p"n^ ~ 8A3/317(2Ai - n^) 

-8A3/3"r2(2Ai - _ l6/3/3"2f73(2Ai - 0^) + 4/3"17(2Ai - r22)2 



SAlVt{l - 13' f + IQAlP^Vt^ - I6A3/3/3" 172(2^1 - 172-) _^ 4/3"2f]2^2Ai - l]^)^ 



= 0. 



(40) 



The real and imaginary parts of the relative transverse shear modes are shown in the Fig. 1 for several values of 
and 77' = = j3" (here and below we scale 77 in units of psi^s)- The real part of the low frequency rotational modes 
is diminished as rj is increased; the high frequency modes are unaffect except in the stong coupling limit 77 > 50 [the 
strong and weak coupling limits are discussed after Eq. (|l^)]. The damping of the modes is maximal for 77 = 1 and 
decreases to zero for 77 —s- and r/ — > 00. Note that in the limiting cases the vortices are locked either in the superfluid 
(77 — > 0) or the normal component (77 00) and hence the damping is ineffective. The transverse shear modes are 
stable for arbitrary values of the eccentricity of the spheroid. 



2. Relative toroidal modes 



These modes correspond to |m| = 2 and the motions in this case are confined to the planes parallel to the equatorial 
plane. The components of Eq. (Cq), which are even in index 3, are: 



^This is easy to prove directly from Eq. (pq). Write the dispersion relation as f{y?) = 0. Then show that (i) f{\^) ±00 
as A^ ±00, (ii) /(O) > 0, and (iii) the two extrema of /(A^) are both at A^ < 0. Thus, the zeros of /(A^) are all at A^ < 0, 
so a = iA must be real. 
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A'f/3;3 


= i^pG)- 


-Mn 


-2A3 


t/33-2f];3"A[/3;3, 






- 2n\U2,i 


- (ttGp)- 






- 2Ai)Un - 2f7/3AC/i;i - 


2r2/3'At/2;i 




+ 2n\Ui,2 


= (^Gp)- 






- 2Ai)U22 - 2npXU2;2 -f 


2n(3'XUi.2 




- 2nxu2,2 


= (-2Ai 




)f/l2- 


^ 2n(3XUi,2 - 2n(3'XU2:2, 






+ 2nxui.i 


= (-2Ai 




)C/21- 


-217/3AC/2;i + 2f7/3'AC/i;i. 





(41) 
(42) 
(43) 
(44) 
(45) 



We add Eqs. (^^ and (|4^) and subtract Eqs. (^2|) and (43) to find the following coupled equations for the toroidal 
modes (note that Ai = A2 for Maclaurin spheroids) 

(a^ + 2n/3X + AAi - 217^) {Uii - U22) - 4f^A(l - j3')Ui2 = 0, (46) 
A^ + 2^i3X + AAi - 29?) U12 + nX{l - /3')(C/ii ~ U22) = 0. (47) 



The characteristic equation for the toroidal modes is 

A'' + 4pX^n + X^iSAi + Ai3^n^ - 8i3'n^ + 4/3'2f]2) 

+X{l6Aij3n - Spn^) + IQAj - leAiil^ + An^ = O. (48) 

In the frictionless limit this can be written as 

(A^ + 4Ai - 217^)2 + 4f72A^ = 0, (49) 

which is factorized by writing A = ia. The two solutions are then 



(71,2 = ± v/4Ai - n^, (50) 

and there are two complementary modes which are found by substituting — il for Q. The modes are always real 
because AAi > H.^ for incompressible Maclaurin spheroids. This result is in contrast to the modes of ordinary 
Maclaurin spheroids which become dynamically unstable at 4i?i2 = fi^, e = 0.953. Our model of the superfluid 
Maclaurin spheroid also becomes dynamically unstable at the same point, but only via the toroidal modes derived 
from the perturbation equations for Vij^ just as in EFE. 

The real and imaginary parts of the dissipative toroidal modes are shown in the Fig. 2, for the same values of rj as in 
Fig. 1. The real parts of the modes tend towards each other and merge in the large friction limit. Note that there are 
no neutral points associated with these modes and the necessary condition for a point of bifurcation is not satisfied. 
The damping of the modes is maximal, as in the case of the transverse shear modes for 77 = 1, and decreases in both 
limits of ry ^ and ry — > 00. In contrast to ordinary Maclaurin spheroids, which become secularly unstable at the 
bifurcation point where 2_Bi2 = and e = 0.813, the new toroidal modes are stable at all values of the eccentricity. 
Our main conclusion is that the toroidal modes associated with the relative motions of the superfluid and the normal 
components always remain stable for incompressible Maclaurin spheroids. In the case of the compressible Maclaurin 
spheroids the point of the onset of secular instability may vary as a function of the adiabatic index (in the case of 
a polytropic type of equation of state) and hence the conclusions reached above should be verified for these models 
separately. 

3. The relative pulsation mode 

To find the pulsation modes, which correspond to m = 0, we first add Eqs. (p2|)-(p3|) and subtract from the result 
the Eq. (pl|). In this manner we find that 



^aV2 + npx + 2Aij {Uii + U22) 

+2nx{i - p'){Ui,2 - U2;i) - (A^ + 4A3 + 2np"x)U33 = 0. (51) 

Subtracting Eqs. ( |45| ) and (Q one finds 

(a^ + 2n^x) {Ui,2 - U2;i) - nX{l - P'){Uii + U22) = 0. (52) 
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Equations (pl\) and (p2) can be further combined to a single equation: 



^A^ + 2f]/3A - 2n^ + AAij (A^ + 2f7/3A) + AQ^X'^il - /3')^J {Un + U22) 
-2 + 2n(3X] (x^ + 2nXf3" + AAs]] U33 = 0. 



The solution is found by supplementing these equations by the divergence free condition 

^11 U22 U33 _ ^ 
af 

or, in terms of the eccentricity e = 1 — a|/a^, 

(f/ll + C/22)(l-e') + t/33 = 0. 



The third order characteristic equation is 
(3-2e2)A3 



8A3 — 0^3 

5//,2o2 



{8f3n + 4(3"n - 4(3e^n - 4:P"e^n)X^ + (AAi 
-2f72 + Ap^n^ - 8f3'n^ + Ap'^n^ + 8(3(3"^^ - 8l3p"e^n^)X 
-8AiPn + 16^3^17 - leAsPe^n - 4/3rj3 = 0, 



where the trivial mode A = is neglected. In the frictionless limit we find (A = ia as before) 



a = ± 



2rj2 + 4A^ +8^3(1- £2 



(3-2e2) 



1/2 



(53) 



(54) 



(55) 



(56) 



(57) 



The pulsation modes for a sphere follow in the limit e, — 0: for a sphere Ai/iirpG) — 2/3, and Eq. ( ^7| ) reduces to 
CT^ ~ 8/3 [a is given in units of (TrpG)^/^]. This result could be compared with the pulsation modes of an ordinary 
sphere: ct^ = 16/15. Thus a superfluid sphere, apart form the ordinary pulsations, shows pulsations at frequencies 
roughly twice as large as the ordinary ones. 

The real and imaginary parts of the dissipative pulsation modes of a superfluid Maclaurin spheroid are shown in the 
Fig. 3. The real parts of the modes are weakly affected by mutual friction and closely resemble those of an ordinary 
Maclaurin spheroid in the frictionless limit. These are located, however, at higher frequencies. The symmetry of the 
damping rate as a function of 77 observed for the transverse shear and toroidal modes is again observed in Fig. 3. 
Note that the results in Fig. 6 were obtained in the case (3" = 0. The pulsation modes of the superfluid Maclaurin 
spheroid are stable, as is the case for the ordinary Maclaurin spheroids. 



B. Modes of superfluid Jacobi eUipsoid 



The sequence of the Jacobi ellipsoids emerges from the Maclaurin sequence at the bifurcation point e = 0.813 
via a spontaneous breaking of symmetry in the plane perpendicular to the rotation (oi 7^ 02). The superfluid 
equilibrium figures are again identical to their ordinary fluid counterparts and the defining relations a\a\Ai2 = 03^3 
and = 2i?i2 are unchanged. Ordinary Jacobi ellipsoids are known to be stable against second order harmonic 
perturbations while they become dynamically unstable against transformation into Poincare's pear shaped figures 
through a mode belonging to third order harmonic perturbations. If the sequence of Jacobi ellipsoids is parametrized 
in terms of the variable cos~^ (03/01), it is stable between the point of bifurcation from the Maclaurin sequence, 
cos^^(a3/ai) = 54.36, and the point where Poincare's figures bifurcate, cos^^(a3/ai) = 69.82. Here, by an explicit 
calculation, we verify that superfluid Jacobi ellipsoids do not develop new instabilities via second order harmonic 
modes of the relative displacements. 



1. Relative odd modes 



The treatment of the oscillations of the Maclaurin spheroid of the previous sections can be readily extended to the 
Jacobi ellipsoids by lifting the degeneracy in indexes 1 and 2 and imposing Ai ^ A2. The equations odd in index 3 
are 
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A'C/3;2 = -2A3C/32 - 2nf3"XU3,2, 
X'^Ui,3 - 2nXU2,3 = {-2Ai + n^) - 2Q.pXUi,-i - 2f7/3'AC/2;3, 
A2c/2;3 + 2nXUi.3 = {-2A2 + n^) c/23 - 2VLpXU2;-i + 217/3'AC/i;3. 
Combining Eqs. ( ^8|) and ( |60| ) and, similarly, Eqs. ( |59| ) and (61), after some manipulation we find 



A^ + 2^2/3" a) (a^ + 2r2/3A) + 2 (^A^ + 2r2/3A) A3 + {x^ + 217/3"a) (2Ai - n^) 
-2nX{l - 13') (x^ + 2r2/3"A + 2A3) C/23 = 0, 



C/' 
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{>? + 2^2/3" a) (a^ + 2n(3X^ + 2 (a^ + 2r2/3A) A3 + [x^ + 2f7/3"A^ (2A2 - O^) 



C/23 



-2r2A(l - /?') (a^ + 2r2/3"A + 2A; 



0. 



Equations (|62|) and ( |63| ) are sufficient to determine the symmetric parts of the virials, and any two of Eqs. 
can be used to find the antisymmetric parts. The sixth order characteristic equation is 



(58) 
(59) 
(60) 
(61) 



(62) 
(63) 



A* 



A^ + 4f7(/3 + p")X^ + 4A3 + 2f7(l - p'Y + A^V{P + I3"Y - {2A^ - n^) - {2A2 - n^) 

+ \l6A3/3n + 8A3(3"n + 8P"n^{l - p'f + + 

-2/317(2^1 - n^) _ 2(3^{2A2 - - 4/3"17(2Ai - if) - 4/3" 1^(2^2 - ^'^] 

+ UaI + 8A3r2(l - /3')^ + 16A3/3^r22 _^ 32A3/3/3"f72 + 8/3"2rj3(l - /3')2 + l6/3^/3"2rj4 

-2^3(2^1 - O^) - 8/3/3" r22(2Ai - _ 4,j3"'^n'^{2Ai ~ VY) - 2^3(2^2 - 
-8/3^"^^2(2A2 - 1^2) - 4/3"2f]2(2A2 - ^ (2Ai - r22)(2A2 - r^^)] A^ 

+ [l6A^/317 + 16A3/3"1}2(1 - /3')^ + 32A3/32/3"173 - 4A3;3f2(2Ai - 17^) - 4^3/31^(2^2 - Vi^) 
-4A3/3"17(2Ai - _ 8/3/3"2rj3(2Ai - O^) - 4A3/3"1)(2A2 - 17^) 

A 



= 0. 



n") 



"^Vt^{2A2 - Vt^) + Al3"Vt{2Ai - 9?){2A2 - Vi^ 
%Al^{\ - /3')^ + l&Ali3^^^ - 8A3/3/3" 17^(2^1 - Vi^) 

-8A3/3/3"rj2(2A2 - 1)2) + 4/3"2r22(2yli - n^){2A2 - r?^) 



(64) 



The real and imaginary parts of the dissipative odd parity modes are shown in the Fig. 4. The Jacobi sequence 
is parametrized in terms of cos~^(a3/ai) starting off from the point of bifurcation of the Jacobi ellipsoid from the 
Maclaurin sequence. The low frequency mode resembles the rotational mode of the ellipsoid; its frequency decreases 
with increasing friction. One of the remaining two distinct high frequency modes is almost unaffected by the dissipa- 
tion, while the other is suppressed close to the bifurcation point in a monotonic manner. The damping rates of the 
odd modes are maximal at 77 = 1 and tend to zero for both large and small friction. The modes are damped along 
the entire sequence; hence, we conclude that superfluid Jacobi ellipsoids are stable against the odd second harmonic 
modes of oscillations. 



2. Relative even modes 



The explicit form of the even parity modes for the Jacobian sequence is 

A'C/3;3 = {TTpG)-^5tV - 2A3U33 - 2np"XU3-3: (65) 

A^C/i;! - 2r2AC/2;i = {7rGp)-^dil - 2AiUii + Q^Un - 2n(3XUia - 2rj/3'AC/2;i, (66) 

X^U2;2 + 217 AC/1.2 = {TrGp)-^Sn - 2A2U22 + ^^U22 - 2ni3XU2;2 + 217/3' AC/i;2, (67) 

A^C/i;2 - 2nXU2,2 = {^^ ~ 2Ai)Ui2 - 217/3AC/i;2 - 217/3' AC/2;2, (68) 

A^C/a;! + 217AC/i;i = (17^ - 2A2)Ui2 - 217/3AC/2;i + 217/3'AC/i;i. (69) 
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These equations can be reduced to a simpler set of equations through manipulations which eliminate the variations of 
the pressure tensor. Explicitly, in the first step we subtract the Eqs. (|6^) and (|67|); in the second we sum Eqs. (|6^ ) 
and (pTI) and subtract twice Eq. (psj). The result is 



(aV2 + npx -n^ + 2Ai)Uii - + nf3x -n^ + 2A2)U22 

~2r!A(l - ;3')C^i2 = 0, 

(aV2 + npx -n^ + 2Ai)Uii + {x^/2 + npx -n^ + 2A2)U22 

-(A^ + 2Vip"X + 4A3)C/33 + 2f^A(l - P'){Ui,2 - U2a) = 0. 



Further we add and subtract Eqs. ( |68| ) and (|6^) to find 



A^ + npX ~ 4Bi2 + 2{Ai + A2) U12 + n{l ~ p')X{Un - U22) = 0, 



(A^ + 2n[3X){Ui,2 - U2;i) - n{l - P')XiUn + U22) + 2(Ai - A2)U-_ 
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0. 



(70) 
(71) 

(72) 
(73) 



Equations ([70|)-(|73|), supplemented by the divergence free condition, Eq. ([m]), are sufficient to determine the modes. 
The characteristic equation is of seventh order, excluding the trivial root A = 0; in the frictionless limit the charac- 
teristic equation is of third order. The explicit form of these equations is cumbersome and will not be given here. 

The real and imaginary parts of the dissipative even parity modes are shown in Fig. 5. For each member of the 
sequence, the eigenvalues of the two high frequency modes are suppressed and that of the low-frequency mode is 
amplified as the dissipation increases. As in the case of the odd modes the damping rates of the even parity modes are 
maximal at 77 = 1 and tend to zero for both large and small friction. The damping rates are again positive along the 
entire sequence and we conclude that superfluid Jacobi ellipsoids are stable against the even parity second harmonic 
modes of oscillations. 



C. Modes of superfluid Roche ellipsoid 

In this section we extend the previous discussion of isolated ellipsoids to binary star systems, and consider the 
simplest case - the Roche problem. The classical Roche binary consists of a finite size ellipsoid (primary of mass M) 
and a point mass (secondary of mass M') rotating about their common center of mass with an angular velocity fi. 
The new ingredient in the problem of the equilibrium and stability of the primary is the tidal Newtonian gravitational 
field of the secondary. Place the center of the coordinate system at the center of mass of the primary with the a;i-axis 
pointing to the center of mass of the secondary and xa-axis along the vector Jl. The equation of motion for a fluid 
element of the primary in the frame rotating with angular velocity Q. is, then (FEE, Chap. 8, Sec. 55) 

dpc, d{dp + ^') 1 a|j^xx|2 

PaUaUa,i = — Pa 7^ ^ -zRa ^ h 2pa€ilmUa,l'^im + i'af3,i, (74) 

where the tidal potential of the secondary, up to quadratic terms in Xi / R, is 



GM' ( xi 2x^ 



(75) 



R \ R 2i?2 

The modified Keplerian rotation frequency for circular orbits, which is consistent with the first order virial equations. 



IS 



r!^ = (l + P)0o(l + 5), (76) 

where P — M/M' is the mass ratio, 0o — GM' /R^ is the tidal potential at the origin of the primary, and 5 is the 
quadrupole part of the tidal field. The latter correction to the Keplerian frequency is maximal at the Roche limit 
where (5 ~ 0.13 [Q. For the sake of simplicity this correction is dropped in the following, as it does not enter into the 
analysis of the stability of the Roche ellipsoid for displacements that are linear functions of the coordinates. However, 
the relation between the frequency f2 and the orbital separation is now determined within an accuracy 5 <^1. As in 
the case of the solitary ellipsoids, we find that the equilibrium figure of the superfluid Roche ellipsoid is identical to 
its ordinary fluid counterpart. 

To treat Roche ellipsoids we modify Eq. (^) to 
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where all frequencies are measured in units of {npGY^^ . Equation 
modes of oscillations of Roche ellipsoids. 



+ [Q^ - MU^, - n^6aU3j + 3cj)o6aUij, (77) 
) is appropriate for finding the second harmonic 



1. Relative odd modes 



The equations determining the modes even and odd in index 3 form separate sets. We start with the modes 
belonging to / = 2 and to = — 1, 1 displacements, which are odd in index 3; for these, 



A^f/g.i = -(2^3 + MU31 - 2r!/3"AC/3;i, 

A2[/3;2 = -2(^3 + MU32 - 217/3" AC/3;2, 



A^C/i;: 



+ 20o) Ui3 - 2rj/3AL/i;3 - 2n/3'XU2-. 



A'C/2.3 



217A[/2;3 = {-2Ai 

2VLXUi,3 - (-2A2 + 0o) U23 - 2npxu2,3 + 2np'\Ui.,3. 

On combining Eqs. (78) and ( pO| ) and, similarly, Eqs. ( |79| ) and ( |8l] ) we obtain 

(A^ + 2f7/3"A + 2A3 + 0o)(A2 + 217/3A) + (A^ + 2Q.p"\){2Ai -9?- 20o) 



(78) 
(79) 
(80) 
(81) 



13 



-2Q.{l-(i')\ 



A^ + 2n[3"\ + 2A3 



U32 = 0, 



(A=^ + 2f7/3"A + 2yl3 + 0o)(A^ + 2f7/3A) + (A^ + 217/3" A) (2^2 - + 0o) 



23 



217(1 -/3')A A^ + 217/3"A + 2A3 
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0. 



(82) 



(83) 



The Uij are symmetric under interchange of their indexes, and Eqs. ( p2D and (83) completely determine the modes 
[any two of Eqs. |78|-|81| may be used to find the antisymmetric parts of t/i^]. The real and imaginary parts of the 
dissipative odd parity modes are shown in Fig. 6, for the case of an equal mass binary {P = 1). The relative modes 
of Roche ellipsoids for other values of the mass ratio display behavior similar to the P = 1 case. The Roche sequence 
is parametrized in terms of cos~^(a3/ai). The low frequency mode resembles the rotational mode of the ellipsoid; its 
frequency decreases with increasing friction. The high frequency modes tend towards each other and merge in the 
limit of slow rotation; in the opposite limit the modes remain unaffected by the dissipation. There are three distinct 
rates for the damping of oscillations. These are maximal at 77 = 1 and tend to zero in both limits of large and small 
friction, as was the case for the Maclaurin and Jacobi ellipsoids. The damping rates are positive along the entire 
sequence. We conclude that superfluid Roche ellipsoids do not develop instabilities via the second order harmonic 
odd modes of relative oscillation. 



2. Relative Even Modes 



As is well known, Roche ellipsoids develop a dynamical instability via the second order even parity modes beyond 
the Roche limit, which is the point of closest approach of the primary to the secondary. If viscous dissipation is allowed 
for, Roche ellipsoids become secularly unstable at the Roche limit via an even parity mode and before dynamical 
instability sets in. We have seen that Maclaurin spheroids do not develop any instabilities (i.e. neither dynamical 
nor secular) via the modes associated with the Uij in the presence of superfluid dissipation. The extension of the 
theory above to superfluid Roche ellipsoids, as we show now, does not reveal any new instabilities in the presence of 



superfluid dissipation, again in contrast to the analysis based on the ordinary viscous dissipation. 
The explicit form of the even parity modes for the Roche sequence is 

A'C/3;3 = {^pG)-^5tl - (2A3 + 0o)f/33 - 217/3" A{/3:3, (84) 

A^[/i;i - 217AC/2;i = (7rGp)"Mn - (2Ai - 17^ - 2(/)o)C/ii - 217/3AC/i.i - 217/3' AC/2;i, (85) 

\^U2,2 + 2n\Ux,2 = {-KCpy^Stl - {2A2 + MU22 + ^^U22 - 2npXU2;2 + 217/3'AC/i;2, (86) 

X^Ui.2 ~ 217AC/2;2 = (17^ + 2fi~ 2Ai)Ui2 - 217/3AJ7i;2 - 217/3'At/2;2, (87) 

A^[/2;i + 217AC/i;i = (17^ - 00 - 2A2)t/i2 - 217/3A[/2;i + 217/3' AC/i;i. (88) 
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These equations can be reduced to a simpler set of equations through manipulations which ehminate variations of 
the pressure tensor. Using the symmetry properties of the virials we first subtract Eqs. (|85| ) and (|86|), then sum Eqs. 
( ^ ) and (|8^) and subtract twice Eq. (^) to obtain 

(AV2 + ni3X ~ 200 + 2Ai)Uii - (AV2 + n(3X -n^ + (l3o + 2A2)U22 

-2r!A(l - /3')C/i2 - 0, 
(AV2 + nf3X -Q^ - 2(j)o + 2Ai)Uii + {X^/2 + r2/3A - f^^ + (/)o + 2A2)U: 



22 



Further we add and subtract Eqs. (|87|) and (pq) to find 



-(A2 + 2r!/3"A + 2(j)o + 4^3)f/33 + 2r!A(l - /3')(f/i;2 - C/2;i) = 0. 



Ui2 + m ~ (3')XiUn - U22) = 0, 



A^ + 2r2/3A + 2(Ai + A2 



2n^ 



(A2 + 2npX){Ui.,2 - U2-i) - n{l - P')X{Uii + C/22) - [30o - 2(yli - A2)]C/i2 = 0. 



(89) 
(90) 

(91) 
(92) 

Equations (|9l|)-(|9^), supplemented by the divergence free condition, Eq. (|5^), are sufficient to determine the unknown 
virials. 

The real and imaginary parts of the dissipative even parity modes are shown in Fig. 7. For each member of the 
sequence the eigenvalues of the two high frequency modes are suppressed and that of the low frequency one is amplified 
with increasing dissipation. In effect these modes merge in the slow rotation limit. The modes do not become neutral 
at any point along the frictionless sequence and, hence, the necessary condition for the onset of dynamical instability 
is not achieved. Note that our model for the Roche ellipsoids is dynamically unstable as is its classical counterpart, 
but via the modes governed by the Eq. ( |l8|) modified appropriately to include the external tidal potential. We do not 
repeat the mode analysis for the virials Vij as it is a complete analogue of the analysis in EFE. As in the case of odd 
modes the damping rates of even parity modes are maximal at = 1 and tend to zero in both limits of large and small 
friction. The damping rates are positive along the entire sequence and we conclude that superfluid Roche ellipsoids 
are secularly stable against the even parity second harmonic modes of oscillations associated with the relative motions 
between the superfluid and normal components. 



IV. VISCOSITY AND GRAVITATIONAL RADIATION 



Above, we neglected viscous dissipation in computing normal modes of a normal fluid-superfluid mixture. However, 
for a single fluid, viscous dissipation is important for understanding stability, for it is responsible for secular instability. 
Although viscous terms spoil the calculation of the modes of uniform ellipsoids from moment equations formally, when 
the dissipative time-scale is long, one can include them perturbatively (e.g. EFE, Chap. 5, §37b). 

The inclusion of viscosity is slightly more complicated for a mixture of superfluid and normal fluid because viscous 
dissipation only operates on the normal fluid. The separation of Eq. (|l^) for the fluid displacements into Eqs. ( p^ ) 
and ( p^ for the common and differential fluid displacements, ^_|_, is possible because the only form of dissipation, the 
mutual friction force, included in Eq. ( pT| ) only depends on d^_/dt. Viscous dissipation depends on only, and, in 
a formal sense, the dynamics no longer separate into the independent dynamics of 

If we assume that the time-scale associated with viscous dissipation is relatively long, then we can include it 
perturbatively. The calculation is a bit more subtle than for a single fluid, because we have to deduce for the 
modes. This brings up an issue that we glossed over earlier, in setting up the calculation of modes in §|l| even when 
computing the modes that arise from the equation for ^_ alone, the equation for must be satisfied, and viceversa. 
The simplest way for this to work is for to vanish when ^_ is nonzero and vice versa. In fact it is easy to show that 
this is a reasonable solution provided that the Eulerian pressure perturbations are 6pa — —£,a,idpa/dxi, a situation that 
arises naturally for adiabatic perturbations, where the Lagrangian pressure perturbations are AaPa = —^ad£a.i/dxi, 
and the perturbations are solenoidal, so that d^a.i/dxi = 0, as is true for all modes considered in this paper. Q To see 
how this works, consider a mode of Eq. (p^), and examine under what conditions Eq. ( p^ will be satisfied. Then, 
using 



^More realistically, one would also need to consider non-adiabatic effects, such as might arise from perturbations from 
/3— equilibrium; see e.g. Lindblom and Mendell These would tend to couple but if small, could be computed pertur- 
batively, as we do here for viscosity, which also couples 
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axi dxi 
and [substituting the definition of and pa = faP into Eq. (p^)] 

^0^-g/^.3,,^(x')^.,(x')A(_2_) (94) 

it is easy to see tliat Eq. ([T^) o nly depends on But since tlie normal modes of Eq. (|6|) have different frequencies 
than normal modes of Eq. (|15|), we must have ^_|_ = when ^_ 7^ 0. Since, in general, 

+ ./jv^- £.N^i+-fs£.-. (95) 

we conclude that, for modes of Eq. (|16|), £0 = /jv£- and — ^fs^-i when Eulerian pressure perturbations are 
given by Eq. (|93|). Similarly, since Eq. ( |l6| ) only depends on for modes with £_j_ 7^ 0, we must have £_ = 
and, therefore, £5 = ^ ^+7 assuming Eq. (p3|). In particular, if the kinematic viscosity ;/ is held constant [see 
EFE, Chap. 5, Sec. 36, Eq. (Ill)], for modes with £^ 7^ 0, the viscous dissipation rate is smaller by a factor of /at 
than it is for a single fluid with the same background and displacement £_|_, as might have been expected c^ualitatively 
(i.e. for small normal fluid density, the viscous dissipation must be diminished). For perturbations with £_ 7^ and 
displacements that are linear functions of the coordinates, we must add 




1 dU,.i 1 dU,. 

5 



to the right hand side of Eq. ( pOD to include the effects of viscous dissipation. 

Perturbations with £_ 7^ emit no gravitational radiation because £+ = for them, and therefore there are no 
perturbations of the quadrupole moment or any other net mass currents associated with them. Gravitational radiation 
is emitted by the modes in which the two fluids move together at the same rate as for a single fluid (e.g. ref. [Q). 
Thus, none of the new modes of a superfluid-normal fluid mixture found here is affected by gravitational radiation at 
all. 



V. CONCLUSIONS 



Despite a number of simplifying assumptions, the study of the oscillation modes of uniform ellipsoids is useful for 
understanding the equilibrium and stability of real neutron stars, at least qualitatively. Moreover, the theory is simple 
enough that it can be extended readily to include modifications to the underlying physics; here, we have considered 
new features that arise because a neutron star contains a mixture of normal fluid and superfluid coupled by mutual 
friction. The theory of uniform ellipsoids is interesting from the viewpoint of mathematical physics because it is 
solvable exactly, and it may also have applications to other physical systems, such as the physics of trapped, rotating 
Bose-Einstein condensates. 

In this paper we have extended previous treatments of the oscillation modes of ellipsoidal figures of equilibrium to 
the case of a mixture of normal fluid and superfluid. The basic equations of motion for the two fluid hydrodynamics 
include mutual friction exactly, since the frictional forces depend linearly on the relative velocity between the two 
fluids, and vanish in the background where the two fluids move together. In addition the fluids are coupled via 
mutual gravitational attraction, which is also treated without further approximations. As a result our relations 
closely resemble Chandrasekhar's tensor virial equations, even though they are intrinsically dissipative due to the 
mutual friction. While we have developed these moment equations for general underlying equilibria, they are most 
useful for perturbations around uniform backgrounds, in which case the moment equations of various orders decouple 
and yield exact solutions for the normal modes. 

Quite generally, there are two classes of modes for small perturbations, one class in which the two fluids move 
together and the other in which there is relative motion between them. The former are identical to the modes found 
for a single fluid. As a result our models of superfluid Maclaurin and Roche ellipsoids undergo dynamical instabilities 
with respect to these modes which are indistinguishable from what is found for their classical counterparts. When 
ordinary viscous dissipation is included they are also subject to secular instabilities related to the modes in which two 
fluids move together. If the kinematic viscosity is held constant, then the rate of viscous dissipation, when computed 
in the "low Reynolds number" approximation, is diminished by a factor /^v, the fraction of the total mass in the 
normal fluid (see however below). 
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The modes involving the relative motion between the fluids are completely new and are shown to be stable along 
the entire sequences of the incompressible Maclaurin, Jacobi and Roche ellipsoids independent of the magnitude of 
the phenomenological mutual friction. These modes also do not become neutral at selected points along any sequence 
and the necessary condition for the point of bifurcation is not achieved. Our main conclusion is that mutual friction 
does not drive secular instabilities in incompressible and irrotational ellipsoids. In addition we find that even though 
the new modes are subject to viscous dissipation (a consequence of viscosity of the normal matter), they do not 
emit gravitational radiation, and are therefore immune to any instabilities associated with gravitational radiation, 
irrespective of their modal frequencies. 

The results summarized above hold within a combined framework of two-fluid superfluid hydrodynamics, Newtonian 
gravity, and the ellipsoidal approximation, as formulated in EFE. Each of these elements of our approach contains a 
number of simplifying approximations which need to be relaxed in realistic applications to neutron stars. For example, 
to treat the mixture of neutron and proton superfluids in the neutron star cores, the one-constituent two- fluid superfluid 
hydrodynamics must be replaced by the hydrodynamics of the multi-constituent superfluid mixtures, in which case 
the mutual entrainment of the superfluids and deviations from /^-equilibrium must be accounted for (see Ref. for 
a discussion of these effects and their impact on neutron star oscillations within the real energy functional method) . 
We anticipate that these effects can be incorporated in the tensor virial approach in a perturbative manner and the 
results of previous sections will hold in leading order of the perturbation expansion. On the other hand, relaxing the 
incompressible approximation and, hence, including the partial pressures of the superfluid and the normal fluid, will 
lead to non-perturbative effects as the pressure terms significantly alter the balance between gravitational attraction 
and centrifugal stretching. As is well known, the points of the onset of the dynamical and secular instabilities of 
compressible ellipsoids depend on the adiabatic index of the underlying polytropic equation of state. As noted in 
Sec. 3, the conclusions reached with respect to the stability of the incompressible ellipsoids should be verified for 
the compressible models anew. The differences in the pressures (or equations of states) of the normal and superfluid 
phases, however, are typically small in neutron stars, since the condensation energy is negligible compared to the 
degeneracy pressures of interacting Fermi-liquids. The coupling between the partial pressures of the noraml fluid and 
superfluid, therefore, can be treated perturbatively. 

Another effect that needs to be included is the density dependence of the mutual friction coefficients and the kine- 
matic viscosity. For example, in neutron stars the kinematic viscosity is density dependent in general, explicitly as a 
result of to the density dependence of the phase space of normal quasiparticles undergoing collisions, implicitly because 
of the density dependence of the in- medium scattering amplitudes (see for further details Ref. |p^ ). The ramification 
for comparison of the secular instabilities of the normal fluid and superfluid ellipsoids is that the modiflcations of the 
rate of the viscous dissipation will depend in a non-trivial manner on the fraction of the normal fluid in the system. 

The entrainment, /3-noncquilibrium, compressibility, e. t. c. will couple the relative and center-of-mass modes in 
general. Such a "mixing", as discussed in Sec. 5, implies ^_|_ 7^ for the relative modes and, similarly, ^_ 7^ for the 
center-of-mass modes. Therefore, the mutual friction might tend to drive the center-of-mass modes secularly unstable; 
if they emit gravitational radiation, the mutual friction will suppress the gravitational radiation induced instabilities. 
One of the important issues to be addressed by the future work is the magnitude of the "mixing" of the modes for 
general equilibria and the corresponding times scales. 
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APPENDIX A: "VIRIAL" EQUATIONS AND PERTURBATIONS 

The two fluids need not occupy the same volume, and we shall suppose that fluid a occupies a volume Va- Taking 
the zeroth moment of Eq. (|l|) amounts to integrating over V^; doing so, we obtainQthe "flrst order 'virial' equation"^ 



^We assume that pa = on the bounding surface of Va. 

^We put the word "virial" in quotes because the equations are dissipative. 
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d 

di 



j d^XpaUa.i 


1 — '^^ilm^jn 1 d 















QiQj) I d^xpaXj 



d'xp^^ + I d'xF^p,,. (Al) 

Apart from inertial forces, which do not couple the two fluids, there are two forces that do couple them: gravity and 
friction. The net, mutual gravitational force between the fluids only vanishes if they (i) occupy the same volume and 
(ii) have densities that are proportional to one another (i.e. ps oc pn)- The mutual friction force is nonzero as long 
as the fluids move relative to one another. Thus we see that, for a two fluid mixture, the zeroth moment of Eq. (|l]) is 
not trivial, as it would be for a single fluid (as in EFE). Note, though, that the center of mass motion of the combined 
system is trivial; i.e., if the center of mass of the combined system starts out at x = with zero velocity, it does not 
movej^ 

Taking the first moment of Eq. (^ results in the second order "virial" equation 



d 
di 



I d XpaXjUa^i j — 26^/771 ( / d XpaXjHaJ \ ^ ^i^k-^a,kj 

JVa. J \JV^ J 

+ (1 - 6o,p)M 



where 



— / d X PaXiXj 



Ua= d XPa 



_ 1 / 3 



2 



G f d^x d^x' pa{x)pa{x'){xi - x'^{xj - x'A 



2 Jy^ |X-X'|3 



d^x'paix.)ppix')xj{xi - a;-) 



Mo,. 

MaP,ij =-G I d^x 

JVa, JVfi 

^aP,i]= j d^XXjFai3,i. (A3) 

When there is only one fluid present, this equation reduces to the results found in Chap. 2 of EFE. There are two 
new terms here: there is a term that arises from the mutual gravitational forces of the two fluids {Aiaf3,ij), and also 
a term from the mutual friction {J-^af3.ij)- 

Consider first the variation of the first order virial equation under the influence of perturbations. Most terms are 
simple to compute, but we must take special care in computing 

d^.p^^^l^^SG I d'xf d'-'P^i^)Ppi^')i-^--d^ (A4) 

^X^ Jv^ Jvf3 |X-X'|3 

In computing the necessary variation, think of x as having a label a and x' as having a label /3. It is then easy to 
find that 

3 / \c / \ ^ f d^x'pp{x'){xi - x'i) 



d-^Xpa^ = -G d''xpa{-X.)^a,l{'X.) 

V„ Jv^ OXl Jv^ 



d^x' Pa{x'){x,- x'^) 

1„ „/i3 (A5) 



There is a slight subtlety that has not been stated explicitly in deriving Eq. (Al). The mutual friction force is nonzero 
only in the overlap volume of the two fluids. This restriction is necessary to derive conservation of total momentum for the 
combined fluids. We might be interested in situations involving long range coupling, such as between the core superfluid and 
crustal normal fluid. For such a coupling to occur, we need to introduce long range fields (e.g. magnetic fields) capable of 
transmitting forces between fluid elements at different points in space. 
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which is manifestly antisymmetric on a ^ /3. Assuming Va = Vp = V and pa = /ap(x) in the background equihbrium, 
we can simplify this to 



d'xpo.jf^ = GUf0 / dMx)[f/3,Kx)-a,z(x)]— / 7 '^ (A6) 

Vq '^■^i J V 1^ ~ ^ I 

Gathering terms, we find that the perturbed first order virial equation is 

^^/q J d^xp£,a,i^ = '2eilm^rn-^(^fa J d^XpS,a,l^ + {fl'^ 6,j - ^.flj) fa J d^Xp^aJ 

d f d^x'p{yi'){x,~ x'i) 



'GUfp d'xp{^)[^fjA^) - ^aM]^J^ 



+5 / d''xFap,^. (A7) 



I Jv 

V 

Although we simplified the final answer by assuming that the fluids occupy identical volumes and have proportional 
densities in the background state, we could not have derived the correct perturbation of the first order virial theorem 
if we had not allowed the volumes to differ. 

For uniform ellipsoids, we can simplify the mutual gravitational term further. First, we recognize that 

f d^x'p{^')ix,-x'^ 90(x) 

so we can write the gravitational term generally as 

Up rf'xp(x) [Cpj (x) - Ui (x)] 1^ . (A9) 
Second, recall that the potential at any interior point of a homogeneous ellipsoid is (Theorem 3 in Chap. 3 of EFE) 

(/.(x)=~7rGp('/-^Afea;?Y (AlO) 



fe=i 



where / is a constant; consequently 



^2TrGpAAi. (All) 



dxidxi 

Thus, the mutual gravitational contribution to the equation of motion for the perturbed center of mass is 



2^Gp'A,fafp I d'x(^/3,,:-Ca,»). (A12) 
Jv 



A sufficient condition for this to vanish is 



/ d'xia = I d'xip, (A13) 
Jv JV 

which is just the statement that the perturbations have the same the center of mass, specialized to the case of uniform 
density. 

In Eq. (A7), we did not compute the variation in the last term. To do this, let us write 

FsN,i = ~Ps^sPij{us.j - unj); (A14) 
in the background state, the two fluids move together (and may even be stationary in the rotating frame) so we have 

6 d'^xFap.^ = -ScPj^ [fs d^xp{^)ujsPr3 {ts.j ~ Cjvj)^ , (A15) 

where Sap = if a = /3, 1 if a = S" and (3 = N , and — 1 if a = and (3 = S, and [see Eq. (^] /3y = 
PSij + (3' eijmVm + (/?" — I3)vii'j- It is clear that the centers of mass of the two fluids remain stationary if the fluid 
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displacements are identical. However, there may be other conditions under which they remain stationary. For example, 
the integrated mutual friction force will be zero as long as 

fs [ d''xp{^)Lusl3^j{^s.j - ^N.j) = 0; (A16) 
Jv 

for a background with uniform density, vorticity and frictional coupling coefficients, this is guaranteed if 

d^'x^s^ I d\iN, (A17) 

JV 

which is less restrictive than the requirement of identical displacements. We found the same condition for the v anishing 
of the integrated, mutual gravitational force for perturbations of uniformly dense ellipsoids. Equation (A17) may be 
true, in fact, for all of the perturbations considered in our paper, since both sides may vanish identically. 

Next, consider variations of the second order virial equation. Most of the terms are varied exactly as for single 
fluids; one exception is 

r /■ , d [ d^x' p{y.')(xi- x'Aix, ~ x'A 

iM.,,. = -G/./,{|<,3.,(x)t(x)- / ix_x,ir 

/■ o d f d^x'p(x'){xi - x'Ax', 1 
+ J^d':,pi^M^,(^)~i^A^)]— J j^_^,j3 (A18) 

where we have specialized to backgrounds with proportional densities and identical bounding volumes. The first term 
in the brackets can be combined with 6Ma,ij and we find 

SMa.ij + (I- Sap)SMap.ij ^ -Gfa / d^ X p{x.)£,a{x) - 



V 



dxi 



/■ , d r d-'x'pix'jixi ~ x,)x'. 
-GUp / d'xp{^)[Ui{^)-^pA^)]— / I ',13 (A19) 

The last equation can be written more compactly in terms of the functions 

_ r d^x'p{K'){x,-x'^)(xj~x'^) dVj _ r d^x'p{-K')x'^{x,-x[) 
= ^ / \7. — zjr?. = / — 3Ts y^^^i 



~ Jv |x-x'|3 dx, - Jy |x-x'|3 

defined in EFE, Chap. 2, Eqs. (14) and (27); using them, we find 

SMo^.^j + (1 - 5^p)5Mc.f3,^J - -fa / d^xp^^j—^ + f^fp / d''xp{^a,i - e/3,i)7r^- (A21) 

Jv (^Xl Jy OXlOXi 

We can rewrite this last result using EFE, Chap. 2 Eq. (28) i.e. 
with this substitution we get 

SMaAj + (1 - Sap)SMaPAj = -fa [ d^Xp^a,l^^ 

Jv (^xi 



dxi dxi dxidxi 



+ fafp / d''xp{Ul-^0,l){^-Sl, 



(A23) 



For the uniform ellipsoids [EFE, Chap. 3, Eqs. (125) and (126)], 



3 3 

= a^xj (^Aj - J2 ^Jkxl^ ^ = ^B,jx,x, + a%, (yA, ~ ^ A.^xl ) , (A24) 
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(ttGp) 



1 d'n 

dxidxi 



^ = -2a'^{SijXiAji + SjiXiAji + SuXjAji) 
= -2a^j[Aii{5ijXi + djiXi) + Sux-jAji], 



dxi 



= 2Bij{SiiXj + SjiXi) - 2af6ijAiixi. 



Using these results and the definitions [EFE, Chap. 2, Eqs. (122), (124), and (125)] 



we find 



(A25) 



(A26) 



5Ma,tj + (1 - Sap)5Map,ij 



nGp 



fa f 2BijVa,ij — a^Sij ^ AilVaM 



1=1 



2Aij{Va,ij — V(3.ij) + Sij ^Aji (]/„,!; — Va,ii) 



1=1 



(A27) 



It is possible to write this more compactly, but the form of Eq. ( A27 ) makes clear which terms depend on the 
differences between the displacements of the two fluids and which do not. 
The other new moment we need is 



5Tal3.,ij = S d^XXjFai3,i. 

Jv„ 



(A28) 



Since the two fluids move together in the unperturbed state, this moment is flrst order in the perturbations at the 
largest. We then find 



5 / d^XXjFaf3,i = -Sapfs / d^XpUsXjPik 



d^S,k di 
dt 



N.k 



dt 



(A29) 



For perturbations of uniform ellipsoids, ujs and p are independent of position in the unperturbed background, and we 
may also assume that /3y is constant; then. 



5 j d^XXjFai3,i = ~Saf3.fs^sPik ^ (^j d^xpXj{^s,k - ^N,k)^ - J d^xpuj{£_s,k - £.N,k) 

For backgrounds in which there are no fluid motions, the last term is absent and 

1 dVa^k-j dVp^k-j 



5 j d XXjFa/3,i = -SalsfsP^sPik 



V dt 



dt 



(A30) 



(A31) 



using the definition in Eq. (A26). 

When there are no fluid motions of the unperturbed star in the rotating frame, the second order virial equations 
are 



dt^ 



2eam^rnfs — + fsVs,!] ~ ^i^kfsys,k] + hj5Iis 



fs7^Gp(2B,,Vs,^J - a^S,, ^ AaVs,i 



1=1 



ajfsfNTTGp 



2A,,{Vs,,j-Vn ,ij) + Sij 



1 = 1 



fs^sPik 



V dt 



dt 



IN- 



'dt^ 



2eiim^mfN + ^'^/nVnaj ~ ^i^kfNVN.kj + 5ij5IiN 
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^ 1=1 ' 



a-ifNfsTrGp 



2Aij{VN.ij - Vsij) + 5,j^Aa[VN,u - Vs,ii) 



1=1 



\ dt dt 
We can replace these equations with a different set by defining 

yi-.j = fsVs.iij + InVna-J, Ui-j = Vs,i-j — VN^i-j. 

In terms of these new quantities we find Eqs. dl^) and (pO[). 



(A32) 
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Fig 1: The real (upper panel) and imaginary (lower panel) parts of the relative transverse-shear modes of a supcrfluid Maclaurin spheroid 
as a function of eccentricity for three values of rj = 0.5, 1, 50. The r) parameter is scaled in units of LUgps- The imaginary parts of the 
modes for r; = 50 are magnified by a factor of 10. The grey lines show the frictionless solutions. To relate the /3-coefficients to the rescaled 
/3-coefEcients we have set fs/fN = 0.2. The results are insensitive to the choice of this ratio. 
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Fig 2: The relative toroidal modes of superfluid Maclaurin spheroids. Conventions are the same as in Fig. 1. 
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Fig 3: The relative pulsation modes of superfluid Maclaurin spheroids. Conventions are the same as in Fig. 1. 
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Fig 5: The relative even-parity modes of superfluid Jacobi ellipsoids as a function of cos ^{a-j/ai). Conventions are the same as in Fig. 
1. 
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Fig 6: The relative odd-parity modes of superfluid Roche eUipsoids as a function of cos ^(as/ai) for P = 1. Conventions are the same as 
in Fig. 1. 
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Fig 7: The relative even-parity modes of superfluid Roche ellipsoids as a function of cos ^(as/ai) for P = 1. Conventions are the same 
as in Fig. 1. 
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